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Abstract 


Calculations are made of the approximate hazard due to peak normal 
accelerations of an airplane flying through a simulated vertical wind 
field associated with a convective frontal system. The calculations are 
based on a hazard metric developed from a systematic application of a 
generic math model to 1 -cosine discrete gusts of various amplitudes and 
gust lengths. The math model simulates the three degree-of-freedom 
longitudinal rigid body motion to vertical gusts with simplified 
expressions for (1) fuselage flexibility, (2) the lag in the downwash from 
the wing to the tail, (3) gradual lift effects, (4) an autopilot, and 
(5) motion of an unrestrained passenger in the rear cabin. Airplane and 
passenger response contours are calculated for a matrix of gust 
amplitudes and gust lengths. The airplane response contours are used to 
develop an approximate hazard metric of peak normal accelerations as a 
function of gust amplitude and gust length. The hazard metric is then 
applied to a two-dimensional simulated vertical wind field of a 
convective frontal system. The variations of the hazard metric with gust 
length and airplane heading are demonstrated. 

Introduction 

The National Aeronautics and Space Administration (NASA) has initiated an aviation safety program 
aimed at reducing the accident rate of transport airplanes. One element of this program concerns the 
injuries and fatalities caused by encounters with turbulence. Although there have been many studies on 
airplane gust loads, these studies have primarily been concerned with the structural integrity of the air- 
planes and the resulting certification issues. However, in most accident scenarios, the airplane structure is 
not damaged, but flight attendants and unrestrained passengers are injured or killed (ref. 1). Often the 
airplane unexpectedly encounters a discrete gust of only a few seconds’ duration while flying in otherwise 
relatively calm air. A study that specifically addresses the transient responses of airplanes to realistic 
turbulence is needed. For example, the critical gust amplitudes and gust lengths that induce accidents 
need to be identified. Once these parameters are known, forward-looking sensors such as radar and lidar 
can be better designed to detect and avoid dangerous gusts. 

A three-step procedure is used to study these factors. First, a simulation math model is developed to 
explore the responses of both the airplane and unrestrained passengers to gusts of various shapes and 
sizes. This model is generic with minimal airplane characteristics. However, the model includes repre- 
sentations of such things as nonlinear aerodynamics of the wing, a simple autopilot, rear fuselage bend- 
ing, the lag in the downwash from the wing to the tail, gradual lift effects on the wing, and inelastic 
collisions of an unrestrained object with the interior of the rear cabin. This model was first described in 
reference 2. Although the model structure is generic, for this study the parameters in the model are 
limited to a single configuration for a medium range jetliner. 

The second step is to develop an algorithm to predict peak normal accelerations for a given vertical 
gust profile. A preliminary version of the algorithm is described in reference 3. The present algorithm is 
believed to be more accurate and easier to implement. It is developed by calculating the transient 
response to a 1 -cosine gust for a matrix of gust amplitudes and gust lengths at a cruise flight condition. 
To make the results more applicable to other airplane configurations and flight conditions, the gust 



amplitudes and gust lengths are nondimensionalized by using airplane response factors. Next, an 
approximate acceleration equation is fitted to the nondimensionalized data. This acceleration equation is 
then used in the hazard metric algorithm to predict the peak normal acceleration for an arbitrary gust 
profile. The algorithm first calculates the approximate incremental gust amplitude for the given gust 
profile for five different gust lengths spanning the estimated range of interest. Then the approximate gust 
amplitude and gust length are substituted into the acceleration equation to predict the peak normal accel- 
eration for each of the five gust lengths. The two peak accelerations (positive and negative) for all five 
gust lengths are defined as the hazard. 

The third and final step is to apply the hazard algorithm and acceleration equation to a simulated wind 
field as an example. The wind field was generated by using NASA’s Terminal Area Simulation System 
(TASS), a large eddy simulation of atmospheric convection (ref. 4). The application of the acceleration 
equation to the simulated wind field produces the two peak accelerations as a function of the assumed 
gust lengths. This analysis gives an indication of the gust lengths that are most likely to be a hazard to 
airplanes. Also, the effect of changing the direction of flight through the same wind field is investigated. 

Nomenclature 

A nondimensional gust amplitude 

A' dimensional gust amplitude, ft/s 

a aft normal force factor at aft passenger cabin (eq. (24), positive upward), ft/s 2 

alt altitude, ft 

a n peak normal acceleration at aft cabin location, g’s 
(a n ) + peak positive acceleration from hazard metric algorithm, g’s 
(a n y peak negative acceleration from hazard metric algorithm, g’s 
z-force factor (eq. (4), positive downward), ft/s 2 
Cf) j total drag coefficient of ith surface, nondimensional 
C Dot zero lift drag coefficient of tail, nondimensional 
C L j lift coefficient of z'th surface, nondimensional 
C [ (l lift curve slope of airplane, per rad 
C Lat lift curve slope of tail, per rad 
c w mean aerodynamic chord of wing, ft 
d fan equivalent fall height in 1-g gravity field (see (eq. 28)), ft 

F s force due to flexible fuselage spring, lb 
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F x j component of force along X-body axis (/ = w, t, T), lb 

F z i component of force along Z-body axis (7 = w, t, ta, T), lb 

/ frequency of flexible fuselage mode, 3 Hz 

g acceleration of gravity, 32.2 ft/s 2 

h position of unrestrained mass above aft cabin floor, ft 

h c collision velocity of unrestrained mass with floor or ceiling, ft/s 

I y moment of inertia about y (pitch) axis, slug-ft 2 

k flexible fuselage spring constant, lb/ft 

kj, altitude hold feedback gain in autopilot, rad/ft 
k 

1 pitch rate gain in autopilot, rad/(rad/s) 
k, efficiency factor of tail, nondimensional 
kfe elevator effectiveness, nondimensional 

k e pitch angle gain in autopilot, rad/rad 

L nondimensional gust length 

L' dimensional gust length, m or ft 

in mass of wing body, slugs 

m e equivalent mass of aeroelastic tail structure, slugs 
nii j least-squares slope for z'th gust length and y'th measurement, (ft/s)/ft 
A« a f t delta normal acceleration of aft passenger cabin, g’s 
An C y delta normal acceleration of center of gravity, g’s 
q pitch rate, rad/s 

q dynamic pressure, psf 

R c ratio of longitudinal cabin position to tail length = 0.8 (eq. (24)), nondimensional 
R s ratio of structural tail deflection to tail length = 0.02 (eq. (12)), nondimensional 
S t horizontal tail area, 542 ft 2 
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S w wing area, 1951 ft 2 
t time, s 

At lag in wing downwash reaching tail, s 

Ug component of gust velocity along X-body axis positive aft, ft/s 

u 0 component of inertial velocity along X-body axis, ft/s 

u w component of true airspeed along X-body axis at wing, ft/s 

V true airspeed, ft/s 

W weight of airplane, lb 

w g vertical gust velocity positive up, ft/s or m/s 

(w'g) j /th measurement of vertical gust velocity, ft/s 

w 0 component of inertial velocity along Z-body axis, ft/s 

w w component of true airspeed along Z-body axis at wing, ft/s 

w x horizontal gust velocity along X-Earth axis, positive aft, ft/s 

w z component of gust velocity along Z-Earth axis, positive up, ft/s 

x position of airplane center of gravity along X-Earth axis (eq. (31)), ft 

Xj x location of point of application of /th force component, ft 

Xj x position of airplane for yth measurement, ft 

x t location of tail aerodynamic center in body axis system, ft 

x w location of wing aerodynamic center in body axis system, ft 

y t ri m vertical displacement of horizontal tail due to aeroelastic effects at trim flight condition, ft 

Zj z location of point of application of /th force component, ft 

a angle of attack, rad or deg 

Ay change in length of spring from trim length, ft 

b ap elevator deflection due to autopilot, rad 

b e elevator deflection, rad or deg 

b ec commanded elevator deflection due to simulated autopilot, rad 
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d£ 

da 


downwash parameter, nondimensional 


ri inclination of thrust axis with respect to x-body axis (positive up), rad 

0 pitch angle, rad/deg 

p density of air, slug/ft 3 

a n standard deviation of normal acceleration, g’s 

o' w standard deviation of vertical gust velocity, ft/s 

Subscripts: 

aft in aft passenger cabin 

b body (axis system) 

e Earth (axis system) 

g gust 

1 /th aerodynamic component or z'th measurement in hazard algorithm 

T thrust 

t (horizontal) tail 

trim at trim flight condition 

ta (horizontal) tail aerodynamic 

w wing 

Operators: 

A change 

first derivative with respect to time 
second derivative with respect to time 
fj) () nonlinear wing drag function 
fl 0 nonlinear wing lift function 
Abbreviations: 
ac aerodynamic center 


c.g. center of gravity 
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NASA National Aeronautics and Space Administration 
TASS Terminal Area Simulation System 

Simulation Description and Results 

Math Model 

The equations of motion are written in a body axis system that is free to translate and pitch in a verti- 
cal plane. The airplane is modeled as a wing-body mass connected by a spring to the effective horizontal 
tail mass. The airplane mass is acted on by the wing-body aerodynamic forces, gravity forces, thrust 
forces, the X aerodynamic force of the horizontal tail, and a Z tail force produced by the connecting 
spring. The Z aerodynamic tail force, the connecting spring force, and gravity act on the effective tail 
mass. A diagram of the model is presented in figure 1 . The equations of motion for this system are given 
below. 

Airplane wing body: 


E F xj ~ m s sin (0) = m (u 0 + qw () ) 

i 

(1) 

E F z,i + m S cos ( 0 ) = ™(w 0 - qu 0 ) 
i 

(2) 

E F xJ z i - E F zJ x i = 7 y<l 

(3) 


i i 


where i = w,t,T. 
Flexible tail: 



m e m e 


-a z -g + x t q 


(4) 


where 


a 


z 


_l 

m 


The X aerodynamic forces for the wing and tail in equations (1) and (3) are derived from the lift and drag 


F x,i = c l s A~ c D,i cos(a f ) + C i>f sin(cq)] 


(5) 
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where i = w,t with the nonlinear, table-lookup aerodynamic characteristics of the wing (see fig. 2): 


F L,w ~ fL^ a w) 

C D,w =/o( a H'^ 

Note that a stall break was not simulated for negative angles of attack because of a lack of reliable 
aerodynamic data in that region. The estimated aerodynamic characteristics of the tail were 

F L,t = C La,t a t 
C D,t = C Do,t + k t C L,t 

The force components for the thrust are 

F x j = Thrust * cos (r|) (8) 

F z j = -Thmst * sin (r|) (9) 

The aerodynamic Z forces for the wing and the elastic tail are 

F z ,i = qSj[-C L i cos(a i)-C D i sin (a ; - )] (10) 




where i = wJa. 

The Z force from the tail ( k 7 , ) that is transmitted to the wing body through the simulated spring is 
not the same as the aerodynamic tail force (F z ta ) given in equation (10). It is, instead, equal and oppo- 
site F s in equation (4): 


F Z,t=- F s= -kiytrim + A -V) (11) 

The spring constant, k, in equation (11) and the effective tail mass, m e , in equation (4) were estimated 
as follows. First, a natural frequency of the elastic tail system was assumed; that is,/ = 3 Hz, which is a 
reasonable value for a medium range jetliner. Then a structural stiffness was calculated by assuming that 
the tail would deflect a given ratio, R s - 0.02, of the tail length if loaded with the entire weight of the 
wing body. That is, 


k = (mg)/(-x t *R s ) 


( 12 ) 


The value R s = 0.02 is consistent with unpublished estimates of the change in horizontal tail incidence 
due to structural deformation for a medium range jetliner. Finally, the effective mass is calculated from 
the spring constant and the assumed natural frequency 


m e 


k 

(2n f) 2 


(13) 
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The angles of attack on the wing and tail are, of course, critical to the proper calculations of the 
aerodynamic forces and moments. The angle of attack at the wing was calculated from the following 
equation: 


a w = tan 


-1 


H' , 


(14) 


in which the velocity components for the wing are a straightforward sum of the inertial components and 
the gust components 


^ Hq ~t" li y 


The primed parameters are the gradual lift effects represented by the approximate formulas (ref. 5): 

w' 0 = w 0 |l -0.36e -r/x ° j 

w g =w g( l - e ~ tlXg ) 


(15) 

(16) 


(17) 


with 


^ c w 

2V 

w 

4V 


(18) 


The velocity components for the tail are more complicated (eqs. (19) and (20)), even though gradual lift 
effects are ignored because of the smaller size of the tail and the shorter time constants due to the shorter 
chord of the tail. For example, the gust hits the tail slightly after it hits the wing; see the u g (t - At) and 

w g (t-At) terms in equations (19) and (20), respectively. Also, the downwash from the wing modifies 

the velocity at the tail; see the |l- — |term in equation (20). Of course, the pitching motion of the 

\ da) ' 

airplane induces a vertical velocity (~x t q in eq. (20)) at the tail that is primarily responsible for the pitch 
damping (ref. 6). Adding the flexible tail degree of freedom adds more terms to these basic effects. The 
vertical motion of the elastic tail induces a vertical velocity component at the tail (Ay in eq. (20)) that is 

responsible for the damping of this structural mode. 


u t = u a + u g ( t - At) + z t * q 


(19) 


de 


da 


de 


w t =w Q - — w a ( t - At) +1 w„{t - At)- x t q + Ay 


da 


( 20 ) 
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where the (t - At) notation indicates the parameter is evaluated at an earlier time (with At = (x w -x t )/V) 
to account for the time difference in effects at the tail and the wing. All other parameters are evaluated at 
the current time t. 

The angle of attack at the tail was calculated from the previously mentioned velocity components plus 
two additional terms: 


-1 w t 

a t = tan — 

V u t I 


+ kfr &e + 


2(A y + ytnm) 


( 21 ) 


The second term on the right-hand side of equation (21) is the effect of the elevator (used from trimming 
the simulated airplane and as the autopilot effector), while the third term is the change in angle of attack 
due to the structural deflection of the rear fuselage. The factor of 2 arises from the assumption that the 
fuselage deflection is proportional to the square of the distance from the center of gravity. 

The autopilot control law contained feedback terms for the pitch rate, pitch attitude, and altitude: 


^ ec kqq + kQ(Q ®trim ) (tilt ^hrim ) ( 22 ) 

This commanded elevator was passed through a 14-rad/s first-order actuator model to produce the 
autopilot output, b ap , which was added to the trimmed elevator position 

. e = ®ap + ^ e,trim (23) 

To simulate the collisions of an unrestrained passenger with the ceiling and floor, an auxiliary calcula- 
tion was made. The first step was to calculate the normal acceleration at an estimated extreme aft cabin 
position 


Titt T; (Rc A/ d 1 //( ■ ) Ay (24) 

where R c =0.8 was the estimated ratio of the aft cabin location to the tail location, x t . The ratio is 
squared in the third term because the structural deflection curve of the rear fuselage was assumed to be 
quadratic. Whenever the unrestrained body (passenger) was not in contact with either the floor or the 
ceiling (0 < h < 4), its acceleration relative to the (accelerating) cabin floor or ceiling was 

h = -a aft (25) 

This acceleration was then integrated to produce the relative velocity and position. Inelastic collisions 
with the floor and ceiling were simulated by setting h = 0 whenever 

h< 0 and (h < 0 or h < 0) (26) 

or 

h> 4 and ( h > 0 or h > 0) (27) 

so that there was no rebound of the passenger from the interior boundaries of the cabin. 
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The 4 (ft) constant in equation (27) was a compromise between (1) the distance between a passenger’s 
head while standing in the aisle and the ceiling and (2) the much longer distance from the ceiling to the 
floor. After the instantaneous velocity, h c , at the time the body collided with either the floor or the 
ceiling was determined, it was transformed into an equivalent height of a fall in a normal 1 -g gravity field 
by using 


dfall = 



(28) 


where V a is an approximate term added to account for the fact that the floor (or ceiling) is either 
accelerating away from or toward the passenger when he impacts the floor (or ceiling). 

V a= + 4 Vl A,1 aft| sign (- A «aft ) (29) 


The constant of 4 (ft/s) is based on the assumption that the body of the passenger deforms 3 in. upon 

impact, 5 = — ft, according to the following equation: 

12 


\V a \- = +p)(32 * | A» aft |) j = + ^16*|A« aft | = +4 


(30) 


Gust inputs 


A 1 -cosine gust profile was used to excite the previous math model. The equation used to produce the 
vertical wind was 


w 0 



1-cos 



(31) 


where x is the x position of the airplane in the Earth axis system. A' is the dimensional gust amplitude, 
and L' is the dimensional gust length. A plot of the 1 -cosine gust profile for typical values of A' and L' 
is shown in figure 3 as a function of time. 

A 1 -cosine gust profile was used because it contains the fundamental characteristics of turbulence that 
excite airplanes — the incremental amplitude and the distance over which the incremental amplitude 
changes. That is, the acceleration of an airplane to a given incremental amplitude gust will increase the 
faster the change in the gust occurs. The 1 -cosine gust profile also approximates a discrete gust rather 
than continuous, random turbulence. Large, discrete gusts are what cause most in-flight injuries (ref. 7). 

Simulated Responses to 1-Cosine Gusts 

The previous equations were translated into a Matlab® Simulink® block diagram for solution. The 
model was trimmed to a cruise flight condition (Mach = 0.8 at 30000 ft) by an iterative routine that 
determined the trim thrust, angle of attack, elevator position, and spring deflection to produce zero accel- 
erations in level flight. The simulated airplane characteristics are presented in table 1 . These characteris- 
tics represent a medium range jetliner that can carry approximately 200 passengers. 

Time histories of the simulated transient response to a large, downward gust starting at approximately 
2.5 s are shown in figure 4. The maximum change in angle of attack was about -7° (for a gust amplitude 
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of -150 ft/s), which did not stall the wing because its initial trim angle of attack was slightly more than 
+ 1°. However, as will be shown later, an upward gust of the same magnitude will stall the wing. The 
normal acceleration at the c.g., A n cg , and the acceleration at the aft passenger location, A« a f t , are practi- 
cally equal because the illustrated gust length of 1 000 ft is so long that the pitching motion, A0, is rela- 
tively slow. As will be shown later, for shorter gust lengths, the acceleration at the aft passenger location 
is considerably higher than that at the center of gravity. 

The change in true airspeed, AH, was much less than the maximum gust velocity of 150 ft/s 
because the latter was nearly perpendicular to the true airspeed of approximately 800 ft/s. The structural 
deflection of the horizontal tail, Ay, was relatively small (less than 0.2 ft), indicating that the simulated 
airplane was relatively stiff. Because the change in the acceleration in the aft cabin dropped below -lg, 
the simulated passenger came off the floor and impacted the 4- ft ceiling at approximately 3.5 s into the 
time history (see the time history of h). The relative velocity at impact was approximately 20 ft/s (see /z); 
however, after the acceleration reversed at about 4.2 s, the simulated passenger fell back down to the floor 
with slightly greater impact velocity than that for the first impact. 

Response Contours 

A matrix of gust lengths and gust amplitudes was run through the simulation program to produce time 
histories. A second program then extracted the maximum values of the accelerations, passenger impact 
velocities, and equivalent fall heights. These data were then transformed into contour plots such as those 
shown in figure 5 for the equivalent fall heights of the unrestrained mass (simulated passenger). The 
maximum contours occur for gust lengths of approximately 500 ft. However, they drop off sharply for 
the smallest gust lengths even though the aft cabin accelerations were very high for short gust lengths, as 
will be shown next. The reason for these dropoffs is that for the shorter gust lengths, the accelerations 
were sustained for such short periods of time that the resulting velocities and displacements of the simu- 
lated passenger were relatively small. In fact, for the shorter gusts the simulated passengers did not 
impact the ceiling at all, but only came off the floor a small distance before impacting with the floor 
again. 

Figure 6 presents the contour plots for the normal acceleration at the c.g. As can be seen in the figure, 
the maximum accelerations are a function of both the gust length and the gust amplitude. For example, at 
a gust amplitude of -50 ft/s, the acceleration varies from -1.0 g’s at a gust length of about 600 ft to only 
-0.5 g’s at a gust length of about 2400 ft. However, at veiy short gust lengths, the accelerations drop off 
again because of the gradual lift effects used in the simulation. Note that for positive gusts greater than 
about 100 ft/s, the accelerations are smaller than the absolute value of accelerations for negative gusts 
because the wing begins to stall, reducing the lift. As seen in figure 2, the lift curve slope begins to 
decrease at about a = 5° and reaches a maximum at about a = 7°. Since the trim angle of attack was over 
1°, the acceleration due to the gust would begin to increase less rapidly whenever the gust-induced angle 
of attack was over 4° (approximately 60 ft/s). As mentioned earlier, nonlinear stall effects were not 
modeled for negative angles of attack because of a lack of reliable aerodynamic data in this region; there- 
fore, the accelerations in the negative direction may be too large for large negative gusts. From a hazard 
standpoint, however, these negative acceleration results would be conservative. 

The contour plots for the normal acceleration at the aft cabin location (fig. 7) are practically identical 
to those at the center of gravity for gust lengths greater than 1000 ft. However, at the shorter wave- 
lengths, the aft cabin accelerations are significantly higher. Part of this increase is due to the pitching 
motion of the simulated airplane, but a significant part is due to the excitation of the structural mode of 
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the tail. The combined effect of the pitching and structural mode is enough to offset most of the reduction 
caused by the gradual lift effects. 


Algorithm Development 

Approximate Acceleration Response Equation 

The contour plots for the normal acceleration at the aft cabin location (fig. 7) were nondimensional- 
ized to make them more generally applicable to other airplane configurations and flight conditions 

2 Weight 

(fig. 8). The gust amplitude was nondimensionalized by dividing by the factor 5 — which repre- 

p svc La 

sents the amplitude of a step gust that would produce a 1-g acceleration. As can be seen in figure 8, the 
y intercepts of the ±l-g contours are approximately ±1, confirming the consistency of the present analysis. 

2 YYl 

The gust length was nondimensionalized by dividing by the factor , which represents the forward 

P SC L a 

distance traveled by the airplane in one time constant of the airplane’s plunge mode. 

To develop a general expression for the hazard metric, a straight-line approximation was fit to the 
contours in figure 8. The straight line that was determined was 

A = 0.96 * a n + 0.70 * a n * L (32) 

where A and L are the nondimensional gust amplitudes and gust lengths, respectively. This equation can 
be solved for the (peak aft cabin) acceleration, a n , 


0.96 + 0.70 *L 


It can be seen from figure 9 that the curve fits the contours better for negative gust amplitudes than for 
positive amplitudes, again because of the stalling of the wing for positive amplitudes. The lack of fit for 
positive amplitudes is in a conservative direction because the curve predicts higher accelerations than the 
accelerations of the contour lines (e.g., +2.0-g contour). 

Hazard Algorithm 

Of course, physical gusts do not generally have 1 -cosine shapes; therefore, it might seem impossible to 
apply this curve fit to real gusts. However, real gusts can be approximated by a series of 1 -cosine gusts 
with different gust amplitudes and gust lengths. In this study a least-squares linear fit was used to calcu- 
late the mnning incremental gust amplitudes for five different gust lengths (see fig. 10). The gust lengths 
were selected so as to span the wavelengths that airplanes are most responsive to -50 m (160 ft) to 800 m 
(2600 ft). To calculate the approximate peak accelerations at the aft cabin location, the incremental gust 
amplitudes were then inserted into equation (33) for each gust length. 

A proposed operational algorithm for a turbulence hazard metric is shown in figure 1 1 . The proposed 
metric actually consists of three numbers: the standard deviation of the acceleration, o n , the positive 

peak acceleration ( a n ) + , and the negative peak acceleration, (a n )~. The metric could be reported 

to pilots as, for example, “Turbulence is 0.2 g’s continuous (a„), with peaks to +0.5 g’s ((n„) + ) and 
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-0.7 g’s ((a n ) ).” Because the peak acceleration, regardless of gust length, may be the critical metric 
(ref. 8), the maximum and minimum accelerations for all five gust lengths combined are reported as 
shown in figure 11. In fact, the peak negative accelerations are probably the most critical because if they 
become more negative than - 1 g, unrestrained passengers will be thrown into the ceiling and then eventu- 
ally back to the floor where they sustain most of their injuries. 

Application to a Simulated Wind Field 

Simulated Wind Field 

Figure 12 contains a contour plot of vertical gust velocities in a convective frontal system as simulated 
by TASS (ref. 4). The gust velocities were calculated on a two-dimensional grid with a 25-m (82-ft) 
spacing in each direction. The plot shows that there are concentrated areas of upward and downward gust 
velocities relatively close to each other. The finer detail of the gust field illustrated in figure 1 3 shows the 
vertical gust velocities on two parallel lines 25 m (82 ft) apart. There are significant differences even this 
close together; thus, it would seem that if this simulated gust field is realistic, two identical airplanes 
flying close together would experience significant differences in the accelerations produced. 

Approximate Acceleration Equation Validation 

Figure 14 shows time histories of a full, six-degree-of- freedom nonlinear simulation of a passenger 
transport airplane flying through the simulated turbulence field. The trajectoiy is approximately along the 
same line as the one in figure 13(a); however, because the simulated airplane was not constrained to a 
given line, the gust velocities encountered were not exactly the same. The simulated peak accelerations 
were approximately ±0.4 g’s from the nominal +1 g. 

The time histories in figure 14 can be compared to the predicted approximate peak accelerations for 
the 100-m (330-ft) gust length shown in figure 15. The predicted approximate peak accelerations are, of 
course, not theoretically the same as the instantaneous accelerations of figure 14. Only the peak values of 
the instantaneous accelerations should correspond to the peak values predicted by the hazard algorithm. 
An operational hazard algorithm (fig. 11) would produce the peaks of all five gust lengths combined, not 
just the 100-m (330-ft) gust length used to produce figure 15. However, as will be shown later, the 100-m 
(330-ft) gust length produced some of the largest peak accelerations. In any case, the predicted incre- 
mental peak accelerations in figure 15 are about ±0.4 g’s, which are in good agreement with the peaks in 
figure 14. 

Acceleration Contours 

Contours of the maximums and minimums for all five gust lengths combined are shown in fig- 
ures 16(a) and (b), respectively. The peak positive acceleration was +1.50 g’s and the peak negative 
acceleration was -1.78 g’s. A comparison of the figures indicates the areas of maximum and minimum 
accelerations are often very close together as a result of the proximity of the up and down gusts previ- 
ously discussed for figure 12. 

Acceleration contours for the five individual gust lengths are shown in figures 17(a) through (e). Al- 
though the numerical values of the contours are different for the different gust lengths, the general pattern 
of the contours is relatively constant. The numerical values of the accelerations are generally larger for 
the smaller gust lengths, indicating that the shorter gust lengths are the most hazardous. 


13 



The acceleration contours in figure 17 are for airplane trajectories traveling from left to right or a 
heading of zero. However, the accelerations at any point are generally a function of the heading the 
airplane flies through the turbulence field. The differences in the acceleration for headings of 30°, 60°, 
and 90° are shown in figure 18 for the 100-m (330-ft) gust length. These figures can be compared to the 
corresponding figure 17(b) which is for a heading of 0°. A close examination of individual points in the 
contour shows substantial differences. For example, the maximum acceleration of 1.42 g’s in figure 17(b) 
occurs at x = 4675 m, and y = 4100 m. The corresponding accelerations at the same point in figure 18 are 
given in the following chart: 


Maximum accelerations at x = 4675 m, y = 4100 m 


Heading, deg 

0 

30 

60 

90 

Acceleration, g’s 

1.42 

0.53 

0.14 

0.37 


The differences in the fine detail mean that airplanes flying through the same point in the turbulence field 
can experience significantly different accelerations depending on the heading. However, as will be shown 
in the next section, the peak accelerations over the entire 256 by 256 points are relatively unaffected by 
the heading. 

Effect of Gust Length and Airplane Heading on Peak Accelerations 

The effect of gust length on the peak accelerations is summarized in figure 19. As mentioned earlier, 
the peak accelerations are generally larger for shorter gust lengths. The peak minimum acceleration 
appears to reach a minimum at a gust length of 1 00 m, but this appearance may be due to the single data 
point at a gust length of 50 m. In any case, if the simulated turbulence field is realistic, this analysis 
would indicate that the largest accelerations are produced by gust lengths of 100 to 50 m (330 to 160 ft) or 
less. This result is supported by actual flight data taken during a NASA research flight for the same class 
of airplane as that used in the analysis. This flight took place in the actual weather system that the present 
simulated wind field was designed to model. The largest “hit” in the flight test produced peak changes in 
the vertical acceleration of +0.81 g’s and -1.35 g’s. The -1.35-g acceleration occurred when there was a 
near step change in the vertical gust velocity of -55 ft/s (17 m/s) over a distance of approximately 64 ft 
(20 m). It is probable that this gust induced a significant structural response. 

Note that since the simulated gust field had a grid spacing of 25 m, the shortest practical gust length in 
this study was 50 m; that is, only 3 points in the grid were used to calculate the incremental gust ampli- 
tudes for the 50-m gust length. Gust lengths of less than 50 m could not be analyzed; therefore, a finer 
grid is needed to establish the value of the gust length that produces the largest accelerations. The peak 
accelerations for gust lengths of less than 50 m need to be determined to see whether the curves in figure 
19 have a local maximum or minimum. However, as discussed earlier in relation to the equivalent fall 
heights of unrestrained passengers shown in figure 5, the gust lengths most hazardous to unrestrained 
passengers may be much higher — on the order of 500 ft. As mentioned earlier, the shorter gust lengths, 
although producing large accelerations, act for such a short period of time that the equivalent fall heights 
are smaller. A hazard metric based on equivalent fall heights would be fundamentally different from the 
present hazard metric based on accelerations. 

The effect of airplane heading on the peak accelerations is summarized in figure 20. It appears that 
variations in the peak accelerations over the entire grid are relatively small and insignificant. However, as 
shown earlier, the differences at a given location can be significant. 
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Concluding Remarks 


A turbulence hazard metric has been developed based on the transient acceleration response of a 
medium-sized jetliner to a 1 -cosine gust at a cruise flight condition. To make the hazard metric more 
applicable to other jetliners and flight conditions, the results were nondimensionalized by using airplane 
response factors. The metric is simple to implement and yet captures the first-order effects of turbulence. 
It is a function of the nondimensional incremental gust amplitudes and lengths. The algorithm produces 
the estimated maximum and minimum accelerations at the aft cabin location. These parameters are 
thought to be a good measure of the hazard to turbulence although other measures such as equivalent fall 
height are also attractive. The metric was applied to a simulated gust field of a convective frontal system. 
A least-squares calculation of the amplitude for five different gust lengths (50 to 800 m) was used on the 
simulated gust field. The analysis showed that the highest accelerations occurred for the subject airplane 
when the gust lengths (for the present simulated gust field) are 100 to 50 m or less. A simulated gust field 
with a grid spacing of less than the current 25 m is needed to determine whether gust lengths shorter than 
50 m produce higher accelerations. The analysis also showed that the estimated acceleration at the same 
point in the simulated turbulence field can be significantly different depending on the airplane’s direction 
of flight. 

The analysis also included a calculation of the motion of an unrestrained mass in the aft cabin to 
simulate the motion of an unbuckled passenger or flight attendant. This analysis showed that a victim’s 
impact against the airplane interior could be equivalent to falling over 10 ft for gust amplitudes of 100 ft/s 
and gust lengths of 500 ft. Thus, from the standpoint of the hazard to unrestrained passengers, the 25-m 
grid spacing of the present simulated gust field is adequate. 
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Table 1. Simulated Airplane Characteristics 


Characteristic 

Value 

Weight, W 

1800001b 

Wing area, S w 

1951 ft 2 

Zero lift drag coefficient of tail, C Do t 

0.01 

Tail lift curve slope, C La t 

1.19/rad 

Wing downwash on tail, — 
da 

0.50 

Pitch moment of inertia, I y 

525 000 slug- ft 2 

Efficiency factor of tail, k t 

0.1058 

Elevator effectiveness, k^ e 

0.5 

x location of wing aerodynamic center, x w 

-1.44 ft 

z location of wing aerodynamic center, z w 

0.0 ft 

x location of tail aerodynamic center, x t 

-66.8 ft 

z location of tail aerodynamic center, z t 

0.0 ft 

x location of thrust, Xj 

0.0 ft 

z location of thrust, zj 

0.0 ft 

Thrust inclination angle, e 

O 

O 

o 

Pitch gain in autopilot, 

2.0 rad/rad 

Pitch rate gain in autopilot, k q 

0.3 rad/(rad/s) 

Altitude gain in autopilot, k h 

0.0001 rad/ft 

Flexible tail spring constant, k 

131 880 lb/ft 

Flexible tail effective mass, m e 

371 slugs 

Thrust, Thrust 

171171b 

Trim angle of attack, a trim 

1.33° 

Trim elevator angle, 8 etrim 

-2.10° 

Trim flexible tail deflection, y trim 

0.0295 ft 

Flexible tail natural frequency, / 

3 Hz 

Flexible tail stiffness factor, R s 

0.02 

Horizontal tail area, S t 

542 ft 2 

Wing lift curve slope (linear range), C La 

7.03/rad 

Total lift curve slope (linear range), C L 

7.625/rad 
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Earth axis 



Figure 1. Axis systems. 



Figure 2. Aerodynamic characteristics of wing. 
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Figure 3. 1-cosine gust profile. A! = 75 ft/s, L' = 1000 ft, V = 796 ft/s. 
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Time, s 
(a) Airplane states. 



Time, s 


(b) Airplane and unrestrained mass states. 

Figure 4. Time histories of airplane and unrestrained mass for downward 1 -cosine gust. A' = -150 ft/s, 
L' = 1000 ft. 
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Figure 5. Contours of equivalent fall heights of unrestrained mass in aft cabin for 1-cosine gusts of various gust 
lengths and amplitudes. 



Figure 6. Contours of peak c.g. acceleration response to 1-cosine gusts of various gust lengths and amplitudes. 
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Figure 7. Contours of peak aft cabin acceleration response to 1-cosine gusts of various gust lengths and amplitudes. 



Figure 8. Contours of peak accelerations in aft cabin as function of nondimensional gust lengths and amplitudes. 
Gust amplitude factor = 34. 15 ft/sec, gust length factor = 844.4 ft. 


21 



Figure 9. Comparison of curve fit used to define hazard metric and previous contours for aft cabin accelerations. 
Dotted lines are curve fit A = 0.96 * a n + 0.70 * a n * L. 
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Figure 10. Illustration of method of calculating incremental gust amplitudes for 50-m, 200-m, and 800-m gust 
lengths using least-squares linear curve fit; 100-m and 400-m gust lengths omitted for clarity. 
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Input Algorithm Hazard 

metric 



Figure 11. Flow diagram for operational hazard metric algorithm. 
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Figure 12. Contours of vertical winds (meter/second) from CFD simulation of convective frontal system (ref. 4). 


23 



(a) Line y = 3200 m. 



(b) Line y = 3225 m. 

Figure 13. Vertical winds along two horizontal lines in figure 12. 
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Time, s 


Figure 14. Real time simulation time history for trajectory approximately along straight line used in figure 13(a), 
y = 3200 m. 



Figure 15. Peak accelerations calculated using curve fit for hazard metric using gust length of 100 m. 
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(a) Maximum acceleration. Peak value = +1.50 g’s. 




(b) Minimum acceleration. Peak value = -1.78 g’s. 

Figure 16. Peak acceleration contours for all five gust lengths combined. 
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Y position, 
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X position, m 


(a) L' = 50m. Peak values = +1.50 g’s and -1.25 g’s. 



X position, m 

(b) L' = 100 m. Peak values = +1.42 g’s and -1.78 g’s. 
Figure 17. Peak acceleration contours for individual gust lengths. 
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(c) L' = 200 m. Peak values = +1.26 g’s and -1.32 g’s. 




(d) L' = 400 m. Peak values = +1.02 g’s and -0.98 g’s. 


Figure 17. Continued. 
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(e) L' = 800 m. Peak values = +0.73 g’s and -0.60 g’s. 
Figure 17. Concluded. 
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1.0 



(a) Heading of 30°. Peak values = +1.27 g’s and -1.41 g’s. 



(b) Heading of 60°. Peak values = +1.29 g’s and -1.30 g’s. 

Figure 18. Peak acceleration contours with L' = 100 m and different headings. 
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(c) Heading of 90°. Peak values = +1.45 g’s and -1.67 g’s. 
Figure 18. Concluded. 
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Figure 19. Peak accelerations as function of gust length for simulated gust field. 
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Figure 20. Peak accelerations as function of heading for gust length of 100 m. 
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